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Abstract 

A random walk problem with particles on discrete double infinite linear 
grids is discussed. The model is based on the work of Montroll and 
others. A probability connected with the problem is given in the form 
of integrals containing modified Bessel functions of the first kind. By 
using several transformations simpler integrals are obtained from which 
for two and three particles asymptotic approximations are derived for 
large values of the parameters. Expressions of the probability for n 
particles are also derived. 

/ returned and saw under the sun, that the race is not to the swift, nor the battle to 
the strong, neither yet bread to the wise, nor yet riches to men of understanding , 
nor yet favour to men of skill; but time and chance happeneth to them all. 
George Orwell, Politics and the English Language, Selected Essays, Penguin 
Books, 1957. (The citation is from Ecclesiastes 9:11.) 
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1 Introduction 



The subject of random motion is one on which an enormous amount of math- 
ematical studies have been made. We mention in this respect the classical 
work of Rayleigh, Smoluchowski, Chandrasekhar, and countless others [1]. 
In this paper, we are interested in the specialization of this general notion 
to random walk on a periodic lattice, where a particle makes random jumps 
between neighbouring sites of this lattice. In this respect we refer in par- 
ticular to the pioneering work by Montroll and his collaborators which has 
provided the inspiration for the present work. 

We shall very briefly indicate the method of Montroll's approach, where 
throughout this paper we shall limit ourselves to random walks on one (or 
more) linear (ID) lattice chains. We shall also suppose that the jump prob- 
abilities of a random walker to the left and to the right are equal, and hence 
equal to p = \. Initially, the time is considered to be discrete, which means 
that we consider the situation of the particle after a discrete number of 
jumps ra, which is equivalent to allowing the particle to jump once in every 
unit of time. Montroll et al. [2], [3] now introduce two quantities which are 
of very great importance. These are 

1. P n (£), the probability that the random walker will be at site £ after 
the nth jump. 

2. f n (£), the probability that the random walker will be at site I after 
the nth jump for the first time. 

Of course, it is assumed that before the first jump (n = 0) the particle is at 
the origin {£ = 0). 

The function P n (£) satisfies the following equation 

P n (£) = \P n -l{£ " 1) + \Pn-x{t + I)- (1.1) 

(If at epoch n — 1, the particle is at either £ — 1 or I + 1, it will have a 
probability p = \ to be at i at epoch n. If it is anywhere else at epoch n — 1, 
its chance of being at £ one jump later, is zero.) This equation also shows 
that the random walk, as described above, is a Markoff process, in that the 
state (£) of the random walker at a given epoch n depends only on that at 
one moment earlier. 

Montroll then introduces a generating function U(£, z) = Y^=o Pn{£)z n . 
This function U(£, z) is then calculated explicitly, from which P n and various 
moments over £ can be calculated. For details we refer to [2] and [3]. We 
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also refer to these papers for the treatment of the first passage times f n (£) 
and the corresponding generating function F(£,z) = J2n!=o fn(£)z n ■ The 
quantity f n {£) is the probability of reaching the site I for the first time at 
the nth jump. 

For the sake of completeness we give the explicit expressions for U(£, z) 
and F(£,z): 



(1.2) 

from which explicit forms of P n {£) and f n (£) follow. 

Montroll et al. [3] also present a method of treating the time as a contin- 
uous variable. Then we introduce as fundamental quantities the following 
probability densities: 

P(£, t) dt : the probability density for the random walker to be at £ 
during interval (i, t + dt). 

F(£, t) dt : the probability density for the random walker to arrive at £ 
during interval (t, t + dt) for the first time. 

(1.3) 

Jumps are now taken to occur at random times h,t2,ts, — This implies the 
introduction of the random variables T\ = t±, Ti = ti — t\, . . .T n = t n — t n -i, 
which have the common density ip(t). For ip(t) we take the exponential 
density tp(t) = ae~ at , where a is the average number of jumps made by the 
random walker per unit of time. From this point on we shall concentrate on 
the first-passage probability density function, that being the one which we 
shall need most in in the following applications. 
We also introduce the probability densities 

Mt)=Ht), M*)= /%(*-r)^ n _i(r)dr, n = l,2,3,.... (1.4) 

Jo 

The function ip n (t) can be interpreted as the probability density that the 
nth jump of the random walker takes place in the time interval (t,t + dt). 
We have 

j.n— 1 

Mt) =a n e~ at (n _ 1)r n = 1,2,3,.... (1.5) 
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It can now easily been understood that [3] 

oo 

F(£,t) = Y,fn(t)Mt)- (1-6) 



n=0 



If we use the tp n (t) given above and the f n (£) that follow from the second 
line of (1.2), we obtain 



oo 



FU t) - 2-V^- 1 V m)n m + 1/2) " {at)£+2n (I 7) 

( ' j " n!(£ + l) n (£ + 2n-l)!- (L?j 

where (a) n denotes Pochhammer's symbol defined by 

(a) = 1, (a)„ = a(o + l)---(a + n-l), n = 1,2,3,.... (1.8) 

Comparing the expansion in (1.7) with that of the modified Bessel function 
of the first kind, see [4, Chapter 9], 

n=U v ' 

and using the duplication formula of the gamma function 

n^) = ^-mr(z+ i 1 ), (i.io) 

an explicit form for F(l , i) is obtained: 

F(£,t) = e- Q ^/,(at). (1.11) 

From [4, Eq. 11.4.13] it follows that for I / 

roc /-oo ,7+ 

/ F(£,t)dt = £ e- at h{at)— = 1. (1.12) 

It is of interest to consider the problem of several simultaneous random 
walkers on a lattice chain, and the behaviour in time of their mutual con- 
figuration. It is as an introduction to this work that we shall consider two, 
three, . . ., independent random walkers on separate lattice chains. We begin 
with two random walkers and consider the situation as shown in Figure 1. 

Remark 1. The integers I a and £b may separately assume negative 
values. However, to avoid the use of absolute value signs, we consider only 
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'a 

B ooooooooooooooooooooooooo 

o 



Figure 1: The random walks for A and B. 

positive values of I a and £ B . But all results hold for negative values when 
we replace these quantities by their absolute values. 

We are interested in finding the probability that particle A arrives at £a 
before particle B arrives at £ B . The solution to this problem is an inter- 
mediate result for the treatment of a ID— grid random walk problem of an 
agglomeration of many particles. 

We know that the probability density for A to arrive for the first time 
at £a in the interval (i, t + dt) is 

p. 

F(e A ,t)dt = -je- at h A {at)dt. (1.13) 
It is now obvious that the desired probability can be written as 

/*oo /*oo 

P(k A <t eB )= dt A F{£ A ,t A ) / dt B F(£ B ,t B ), (1.14) 
Jo Jt A 

where t( A is the time that particle A reaches the site £a for the first time, 
and similar for t£ B . The independence of the walkers is expressed by the 
fact that it is the product of two F— functions which is being integrated. 
Using (1.12) we have 

roo rt\ 

P(t tA <tt B ) = l- dt A F(£ A ,t A ) / dt B F(£ B ,t B ) (1.15) 
Jo Jo 

and interchanging the order of integration in this integral, we derive the 
symmetry properties (which are evident from the random walk problem) 

P(t tA <t tB ) = l-P(t tB <t tA ), hence P(U A = t iB ) = ±, ii£ A = £ B - 

(1.16) 
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If I a = Ib w e can also use integration by parts 

/ dt A F{£,t A ) / dt B F{£,t B ) = 
Jo Jt A 

- / / F(£, t)cLt d / F(£, r)dr =\ / F(£, r) dr 
Jo UtA J \~Jt A J Ljo 



i 

2' 



(1.17) 

In this paper we derive asymptotic expansions of P(tg A < tf B ) given by 



P = 



r^e-^I eA (tA) r^e-^h B (t B ). (1.18) 

JO M Jt A *B 



In (1.18) the scale factor a has been absorbed in t A and t B , because of 

dt dat 

t at ' 

We will give one expansion that holds for large values of £ A and one 
for the case that both parameters £ A and £ B are large. We also give an 
expansion that holds just when the sum £ A + £ B is large. 

2 Transforming the integral 

We study the integral (1.18). We use well-known properties of the modi- 
fied Bessel function to transform the double integral in (1.18) into a single 
integral. 

The inner integral in (1.18) can be modified by evaluating 

S e (t) :=£ e - s h{s)-. (2.1) 
Jt s 

where £ = 1,2,... . We use the integral representation (see [4, Eq. 9.6.19]) 



Us) 



1 

7T JO 



e scosd cos nOdO. 



for integer values of n. Integrating by parts we obtain 
-IJ S ) = - f W e scose sme sm£6de. 

S TT Jo 



(2.2) 



(2.3) 



It follows that 

roo 



POO 1 f7t /"OO 

/ e- ps h{s)— = - Main 9 sin £6 ds e - s{p - cose) 

Jt S TT Jo Jt 



e - P t f n s[n s[n ig 



(2.4) 



£_ f 



,t cos 9 



p — COS ( 



de, 
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which holds for p > 1. It follows that Se(t) of (2.1) can be written as 

Se (t) = — r si * 9sinie e tcose dd= ^!_ r cot y sulWe tco S de (25) 

IT JO 1- COS 9 ^ JO 

Using this relation and interchanging the order of integration in (1.18), we 
obtain 



(2.6) 



de. 



P = iA e-% A (t)S eB (t)- 
Jo t 

JaT^q sui£b [f e -^»»j iA (ti 

vr Jo Uo t 

Invoking again (2.3) we obtain 

1 [* sin^sin^B^ [* , D sin 9 X sin ^ 
p= d6 2 — : - — / d0i t- BTTTi TT- 2 - 7 ) 

7T 2 Jo 1 - COS 9 2 Jo (1 - COS 9\) + (1 - COS # 2 ) 

The — integral can be evaluated; see [5, Eq. 3.613(3)]. Another way is to 
use in (2.6) the Laplace integral 



I 



/ o °° e-^l e (t) j=(p+ VP*^) > I > 0, V > 1, (2.8) 



which follows from [4, 29.3.53] by taking a = 1, b = — 1. This gives 

P = 1. J\ot\e sm£ B 9 (p+yV-l) p = 2-cos#. (2.9) 

3 Asymptotic expansions 

We give three asymptotic expansions: 

• one for large I a-, with £g fixed, or small, 

• one for large I a and £b, with I a ~ £b, 

• one uniform expansion in which one or both parameters may be large. 



3.1 The case i A > t B 

We start from (2.9) in the form 



p=- r /(0) e -^w d 5 (3.1) 

JO 



where 



f{0) = cot\9 sm£ B 9, (f)(9) = In (p+\]p 2 - , p = 2-cos9. (3.2) 
First we observe that 

sin# cos^ 
VP - 1 yl + sin 2 ±0 

Hence, 0(0) is an increasing function on [0, tt] with 

0(0) = 0, 0'(O) = 1, 0'(7r)=O. (3.4) 

It follows that 

P ~ - / f(9)e~ eA(p W d 9, (3.5) 
Jo 

where #o is a fixed number in (0, tt), and the error in this approximation is 
exponentially small when Ia is large. 

Carrying out an integration by parts in the form 

„ -i f e ° f(0) 



tt£ a 

leads to 

e=e 

"~ n£ A 

u— u 

where 



Jo m {d - b) 



P - ^< " ' + -4- / ° /i(0)e"^^ dff, (3.7) 

TTIA JO 



We can repeat this procedure, and compute the integrated terms. The terms 
at 9q can be neglected because they give exponentially small contributions 
compared with the contributions from 9 = 0. Note that we cannot take 
9q = tt, because 4>'(tt) = 0. 



8 




Figure 2: Graphs of P(tt A < t( B ) based on the asymptotic approximation (3.9) 
(shown in red), compared with graphs based on the expansion (3.34) (shown in 
blue) . The red and blue graphs deviate from each other because of the failure of 
the red approximations for large values of is- 



In this way we obtain the asymptotic expansion 



7r£ A 



01.02. 

ao + T A + e A + --- 



where, for k = 0, 1, 2, . . ., 
/fc(0) 



Ofc 



d h{e) 
de pipy 



M0) = f(o). 



(3.9) 



(3.10) 



The coefficients with odd indices are zero. This follows from observing 
that f(6) and 4>'(9) are even functions; see (3.2) and (3.3). Hence, fi{9) of 
(3.8) is odd. By using the recursion in (3.10) it follows that f2k{@) is even, 
and that f2k+i{@) is odd. The first non-zero coefficients are 



Ofj 



2£ B , 02 = !Ml-4), 04 = ^(23- 804 + 124). (3.11 



In Figure 2 we compare the approximations based on (3.9) with values 
obtained by using the expansion in (3.34), which holds when £a + £b is 
large. We see that for smaller values of £b the graphs of the asymptotic 
approximation (3.9) are in agreement with the graphs obtained from the 
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expansion that holds when at least one of the parameters I a ox I b is large. 
The failure of the red approximations is due to the failure of the asymptotic 
approximation (3.9) that has been chosen for this case. 

3.2 The case £a ~ (-b-, both large 

We replace in (2.9) t A by I and lg by i + 5. We know that P = \ if <5 = 0. 
We expand (2.9) for large values of £, keeping 5 fixed. We have 

P = | + Pi + P 2 , 

where 



i jT cot \6 (cos 59 - 1) siuiO (p + - 1 



1 



3- y cot ±0 (cos <J0 - 1) e ltt) [ p + yp 2 ^ 



d0 



d0, 



1 /■* 1 
— / coth 

TT JO 

1 Z" 71 
»- / CC)t± 
TT JO 



9sm59cos£9 (p+\p 2 -l) d6 



9 sin 59 e m [p + 



\Jp 2 ~ 



d9. 



Let 



g (0) = cot \6 (cos 56-1), (f)(6) = -id + In (p + ^p 2 - 1 
Then we integrate by parts in the integral for Py. 

and we obtain an expansion as in (3.9), 



1 , h b 2 
bo + j + ^ + 



where, for k = 0, 1, 2, . 



6 fc = 9 1777^' 9k+i{0) = —- [f7 ^, g (6) = g{6). 



'(0) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 
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It turns out that the coefficients with even indices are zero. To verify this 
we can use a similar argument as for the in (3.10). The first non-zero 
coefficients are 



bi = -\5\ h = \5\ & 5 = A^ 2 (4<5 4 - 205 2 -77). 
In a similar way, let h{6) = cot \Q sin 88. Then 

P2 



(3.18) 



1 



ci c 2 
c + j + ¥ + ... 



(3.19) 



where, for k = 0, 1, 2, . . ., 



(o) 



hk+l{e) = WW)' ho(0) = h{0). (3.20) 



It turns out that the coefficients with even indices are zero. The first non- 
zero coefficients are 



co = 5, c 2 = U(5 2 -l), C4 = --U 2 (12<S 4 + 20<5 2 -77). (3.21 



3.3 The case i A + £ B large 

Because of (see (3.4)) 

we have for large values of £a 

p ^2_ r^ e -e A e de = l avcta jB ^ 

where we used [4, Eq. 29.3.110]. 

Observe that this estimate perfectly reflects the properties of P men- 
tioned in (1.16); also, it is less than unity, as the probability P itself is. 
Moreover, in this estimate large values of Is do not disturb the approxima- 
tion. 

The result (3.23) is obtained by combining the dominant behaviour of 
e -^A0(0) near the origin with the complete form smisO, without expanding 
this function. 

We modify the integration by parts procedure of §3.1, by including the 
(possible large) parameter in the "phase function" <p(6). We can do 
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this by writing smisO = 9e i£s6 '. A complication is the pole of the func- 
tion cot \9, which singularity is removable in combination with the function 
sin IbO. 

To perform the integration by parts procedure we proceed in the follow- 
ing way. In (2.9) we can consider £b as a continuous parameter, and we can 
differentiate with respect to Ib- We also observe that P vanishes with £b- 
We have 



8P =i» 



/■7I 

Jo 



d£ B vr 
We write this in the form 

where 
with 



0cot±0e MBW p + \ p 2 - 1 



d6 



o£b vr 

Q= r f(9)e m d0, 
Jo 



f{9) = \0 cot \9, i>{9) = i£ B 9 - £ A Hp + \jp 2 - 1) 
We integrate by parts, starting with 



(3.24) 

(3.25) 
(3.26) 

(3.27) 



f(0) 



(0) 



where 



f(0) 



MO) 



+ 



r fi(9)e m d9, (3.28) 
Jo 



d f(9) 



d9f{9Y (3 ' 29) 

We repeat this procedure, and compute the integrated terms. Again, the 
terms at 9 = ir can be neglected. 
We obtain 

Q ~ d + di + d 2 + . . . , (3.30) 

where 



dk 



fk(0) 

#(o)' 



fk{9) 



d f k -i(9) 
d9 y/(0) ' 



k = 0,1,2, . 



(3.31) 



and fo(9) = f(9). Again, all coefficients with odd index vanish. This follows 
from 



V>'(0) = %£b - U 



sm( 



1 



Ub - £a 



cos ■ 



1 + sin 2 h 



(3.32) 
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which is an even function and f(6) is also even. Hence, f\(0) in (3.29) is 
odd; and so on. 
We have 

1 _ It A + i£ B 23^ + 129j^g + 2€| 

d °-£ A -U B > d2 -Q(i A -u B y> d4 - 60(^-« B ) 7 

(3.33) 

Considering (3.25), taking the real parts of the coefficients and integrat- 
ing the real parts over the interval [0,£b], we find 

P ~ - (e + e 2 + e 4 + e 6 . . .) , (3.34) 

7T 



where 

The first few are 
e = arctan^, 

p IaIb(1 2 a -1%) 



e 2k = f tB d 2k (£' B )d£' B . (3.35) 
Jo 



e 4 



3(^+4) 3 ' 

I^b^a - l|Kg3^ - 354^| + 23£|) 
60(^+^)6 

^aM^i - £ 2 B )(249£ A - IQ79M A £ 2 B + 40630^^ - 10796^^ + 249£|) 



126(^ + 4)9 

(3.36) 

We see that the shown coefficients e 2 , e±, e§ vanish when £a = £b, and that 

in fact e 2 k{(-A^B) = \^^k,0 ~ ^2k(^B,^A), k = 0,1,2, These properties 

are in agreement with the relations for P in (1.16). Because there is no 
symmetry in (2.9) with respect to £a and £b, they do not follow from the 
construction of the coefficients d 2 k and e 2k . 

When we scale the parameters by putting £ B = \£a, we see that the 
shown coefficients obey the relation 

d 2k = O (£ A 2k ) , (3.37) 

uniformly with respect to A > 0. When we write (3.30) with a remainder, 
that is, 

Q = d + d 2 + . . . + d 2k _ 2 + r f 2k (9)e^ d9, (3.38) 

Jo 
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a straightforward analysis shows that similarly 

f2k(0) = O (r A 2k ) , (3.39) 

uniformly with respect to A > and 6 G [0, ir]. This shows the nature of 
the uniform asymptotic expansion of Q, and, after integrating, the nature 
of the expansion for the probability P. 

By expanding the coefficients et in (3.36) for large £ A with Ia 3> £ B , we 
obtain the coefficients of the non-uniform expansion of §3.1. 

4 Three particles and more 

For three random walkers A, B, C the probability integral reads 
P(t tA < t tB <t tc ) = 

roo roo roo 

/ dt A F{£ A ,t A ) / dt B F{£ B ,t B ) / dtcF(£c,t c ), 
Jo Jt A Jt B 

with the density as in (1.13). That is, 



P(te A < k B < t tc ) 

r°° dt A 



e c r^ w r^-u.(u) r^e->°i lc (tc). 

JO tA Jt A *B Jt B tC 

(4.2) 

It gives the probability that particle A reaches site £ A , before particle B 
reaches £ B , while B reaches site £ B , before particle C reaches £c- 

First we observe that the probability for three particles arriving at the 
same site £, that is, £ A = £ B = £c = £ equals J,. This easily follows from 
(cf. (1.17)) 



poo roo r roo "I 2 

/ dt B F{£,t B ) / dt c F{£,t c ) = \ / F{£,r)dT 

Jt A Jt B IJtA 



(4.3) 



Substituting this in (4.2), performing another integration by parts, and using 
(1.12), gives the value Jy. Using the same method we infer that for n particles 
the probability for all n particles arriving at the same site £ equals 

Repeating the steps used for obtaining (2.7), and replacing all Bessel 
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functions by using (2.3), we easily find for (4.3) 

sin 9 sin let 



1 f n ^sin^s 

P{k A <t lB <t lc ) = - d9—— 

7T 5 Jo 1 — 



X 



cost/ 

(4.4) 

sin a singer f T sinr sin I at 



[* , sincr sin£ B ff / >7r , 

/ da / dr ■ 

Jo 2 — cos — cos a Jo 



3 — cos 9 — cos a — cos r 

Evaluating the r— integral gives 

7T Jo Jo 1 — cos # 2 — cos # — COS <7 



jszijj v5^) n 



(4.5) 

where g = 3 — cos — cos <r. 

From the above analysis it is clear how a similar integral representation 
can be obtained for n random walkers A±, A2, . . . , A n . The probability can 
be written in the form of the n— fold integral 

1 r r A sin 9j siniA 9-j 
P(te Ai <t eA2 <...<t £A J = - d9 1 ... d9 n H 1 , 

(4.6) 

where 

n n 

i J i = E( 1 - cos ^) = 2 E sin2 |^> j = l,2,...,n. (4.7) 

k=j k=j 

Integrating the 9\ integral gives 

- l M ^ sin9j sm£ Aj 0j 

,2 PJ 

(4.8) 

where 

n 

p = 1 +p 2 = n - ^cos 9j. (4.9) 
i=2 

4.1 Asymptotic approximations for three particles 

For large values of £a the main contributions to the integral in (4.5) come 
from the origin a = 0, 9 = 0. To see this we observe that 

q + - 1 = 

3 — cos 9 — cos a + \J {2 — cos 6* — cos a) (4 — cos 6* — cos a) = (4.10) 
1 + V9 2 + a 2 + O (V, 0a, a 2 ) , 
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and that 



(q + y/q* - l) = e-^fo+V?^) „ e ~^VW (4 . n) 



as 0, cr — > 0. We also have 

^sin^ „ crsino" 2<r 2 



1 — cos 9 2 — cos 9 — cos <r 8 2 + a 2 

as #, a — ► 0. 

This motivates us to consider as a first approximation 



(4.12) 



P ~±r d e r d a *^ ^_ uve^^ (4 . 13) 

vr 2 Jo Jo 9 a V 2 + cH 

where we have used (4.5), (4.11) and (4.12). 

Next we use polar coordinates for 9 and a by writing 

9 = rcos(f), a = rsmcf), < r < tt, < (j) < ^ir. (4.14) 

We extend the finite square in the (9, a)— plane to the quarter plane and 
obtain 

P~4 [*"&!> r r dr Sln( ^ r C ° S ^ Sln( ^ r ^ sin 2 e ^ . (4.15) 
Jo jo t cos r sin <p 

The r integral can be found in [5, Eq. (3.947)], that is, 

f°° -or ■ \ ■ i \ dr i, a 2 + (b + c) 2 

/ e sm 6r sin cr — = ± In ^ , 4.16 

7o r « + (o — cr 

and can be proved by differentiation with respect to a. We obtain 

p 1 M*. ,. ^+(^cos0 + ^sin0) 2 

P 2/ tannin — - — — — #, 4.17 

7T Jo «1 + («c cos (f) - £ B sin </>) 2 

which can be written as 

P 7 / tannin- --7 —77 #, 4.18 

tt 2 7 1 + w cos 2<p - v sin 2</> r v ; 



where 

u = o«2 i - ^ , i>2 ' u ~~ 0^2 I 71 I 7T m (4.19) 



i 2 c-i 2 B 2£ c h 
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When v is small we can expand 



v sin 2<j> 



1 + u cos 20 + v sin 20 1 t i +UCOS 2</, 

In — : — — = In - 

1 + u cos 2<p — v sin 20 ] 



2E 



sin 20 

1 + MCOS 2(f) 

1 t; 2n+1 sin 2w+1 20 

2n+l 



(4.20) 



^ 2n + 1 (1 + ucos20) 



which gives 



2 ~ t; 2 " +1 J sin 2n+1 20 ;j ^ , 

P~— V / tan0- ^-rrd0- 4.21 

vr 2 ^ 2n + li (1 + ucos20) 2n + 1 Y y ' 

This expansion is useful when I a is large compared with £b and £c- 

The integral in (4.21) can be written in terms of a Gauss hypergeometric 
function, and the sum can be written as an Appell function. This does not 
give further insight, however. We prefer to give a few further estimates. 
For examining the convergence of the series in (4.21), observe that 

sin 2n+1 20 2 . . 

tan 7 -V-tt < 7 \o — tt ) 4-22 

^(1 + ncos20) 2n + 1 ~ (l-u) 2 ™ +1 ' V ' 

with 

2 e A +e B +e c ' {4 - Z6) 

which is bounded away from 0, unless ic is much larger than I a and £b- 

It follows that expansion (4.21) can be viewed as an asymptotic expan- 
sion for small values of v for the right-hand side in (4.18). 

Of further interest is that when u = 0, that is, £b = f-c we can evaluate 
the right-hand side of (4.21) in terms of elementary functions. In fact we 
obtain by using 

^sin 2 "+ 2 0cos 2 "0# = r(n T t| )r j _ ( " + * } , (4.24) 
o T(2n + 2) 

a Gauss hypergeometric function, that can be written as an elementary 
function: 

P ~^2Fi(|,|;|;^)=iarcsin W , v = J^fpT- ( 4 - 25 ) 
When I a = £b = £c this becomes 



L 



P ~ iarcsin- = -, (4.26) 

7T 2 6 



which is the exact value. 
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5 Discussion and concluding remarks 



We have discussed in this paper a method of considering different simulta- 
neous independent ID— random walks. This work has been motivated by an 
attempt to describe the agglomeration of a number of random walkers on a 
linear chain which will be fixed when they come to occupy nearest neighbour 
positions on the chain. In treating this problem it turns out to be possible to 
effect a transformation of coordinates which makes the evolution equation 
become separable, such that we obtain a product of " one-particle" equations 
which can then be mathematically treated as independent random walkers 
as described in this paper. However, it turns out that this separation is 
possible only when the jump probabilities in both directions are equal. This 
is the reason why we have limited ourselves to equal jump probabilities in 
this work. 

For two particles we have given a complete asymptotic description for 
the case when I a and/or £b are large. For three particles we have also 
given asymptotic results, but a full description becomes a very complicated 
matter. 

Very recently a paper [6] has appeared which treats a related problem 
(with discrete time steps) by a different method, involving stochastic matri- 
ces. 
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